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I.  INTRODUCTION 


The  purpose  of  this  report  is  to  derive  a  set  of  equations  describing  the 
motion  of  artillery  shells  for  convenient  solution  by  computer.  The  effect  of 
mass  unbalance  is  taken  into  account.  The  equations  are  set  up  in  an  inertial 
coordinate  system  and  the  effect  of  the  rotation  of  the  earth  is  included. 

The  form  of  the  ballistic  equations  for  rockets  as  formulated  in  University 
of  Illinois  TAM  Report  No.  166  is  preserved  whenever  possible,  and  the  same 
notation  is  used. 

The  method  of  treating  the  rotation  of  the  earth  is  that  used  in  University 
of  Illinois  TAM  Report  No.  204. 


2. 


11.  GENERAL  CONSIDERATIONS 


1.  Equations  of  Motion 


The  artillery  shell  is  a  rigid  body,  and  the  somewhat  complex  theorems 
of  linear  and  angular  momentum  which  must  be  used  in  rocket  theory,  take 
the  following  elementary  form: 

The  vector  sum  of  all  the  exterior  forces  acting  on  a  system  of  particles 
is  equal  to  the  time  rate  of  change  of  momentum  of  the  system. 

The  vector  sum  of  the  moments  of  the  exterior  forces  acting  on  a  system 
of  particles,  taken  about  the  centre  of  mass  of  the  system,  is  equal  to  the  time 
rate  of  change  of  the  angular  momentum  of  the  system.  In  equation  form: 


“i;  *=  ■  “  t 


=  G  +  F 


(1) 


^  dT  =  ^  (2) 

where  ZFu  I]  r  x  F  represent  respectively  the  sum  of 

^  K  ^  K  k 

the  exterior  forces  acting  on  the  shell  and  the  sum  of  the  moments  of  these  forces 
taken  about  the  mass  centre;  Ttf  ,  and  are  respectively  the  momentum 

of  the  shell,  and  the  angular  momentum  of  the  shell  about  the  mass  centre.  G 
is  the  resultant  aerodynamic  force,  F  the  force  due  to  gravity,  and  N  the 
moment  due  to  the  aerodynamic  forces. 


o 
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2.  Coordinate  Systems  and  Basic  Assumptions 

In  establishing  the  equations  of  motion  of  an  artillery  shell,  we  must  apply 
Newton’s  laws  of  motion  in  an  inertial  coordinate  system.  For  purposes  of  this 
discussion,  we  will  consider  an  inertial  system  to  be  one  which  consists  of  the 
normal,  tangent,  and  binormal  to  the  path  of  the  earth's  cent  re  The  particular 
system  to  be  used  will  be  rectangular  X^,  Y^,  coordinates  chosen  in  such 
a  way  that: 

(a)  The  origin  will  coincide  with  the  mass  centre  of  the  shell  at  launch. 

(b)  The  direction  of  aim  of  the  shell  will  be  in  the  plane  =  0. 

(c)  The  axis  will  pass  through  the  "centre"  of  the  earth. 

In  order  to  obtain  the  location  of  the  shell  in  a  system  fixed  with  respect  to  the 

earth,  we  have  to  apply  a  transformation  between  two  coordinate  systems.  These 

transformations  are  derived  in  the  appendix.  Let  us  briefly  describe  the  nature 

of  the  two  coordinate  systems.  The  coordinates  (X^.,  Y^,  Z^  )  (measured  with 

respect  to  the  earth),  are  defined  in  Eqs.  (A-17)  to  be  functions  of  (X  ,  Y  Z  1 

o  o’  o' 

(the  inertial  coordinates),  and  also  of  (®j^>  ),  the  coordinates  and 

direction  of  the  shell  at  launch  measured  in  a  system  fixed  in  the  earth. 

It  may  further  be  convenient  to  express  the  (Xg,  Y^,,  )  coordinates  in 

terms  of  the  longitude,  latitude,  and  altitude.  In  Eqs.  (A-19)  (p„,  <j>  ,  0  )  are 

E  E  E 

expressed  as  functions  of  (X^,  Y^,  Z^  )  .  Here  and  9^  are  related  to  the 
latitude  and  longitude  as  described  in  section  iii  of  the  appendix. 

The  wind  velocities  which  are  used  in  computing  the  aerodynstmic  forces 
must  be  expressed  in  terms  of  the  (X^,  Y^,  Z^  )  system.  However,  they  will 
be  given  in  terms  of  the  earth  coordinates  (0^,,  (js^,  p^  )  ,  and  are  specified  by  a 
direction  ip^  (azimuth)  plus  vertical  and  horizontal  components  P  and  Q 
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Let  us  now  return  to  the  inertial  coordinate  system  and 

talk  about  the  equations  of  motion  in  detail. 


orientation  of  the  X  Y  Z  system 
o  o  o 

given  in  terms  of  the  angles  0.  ,  4>t  . 

J-# 

ij)j^  measured  relative  to  the  (X*,  Y* 
system,  as  described  in  the  appendix 


X*  =  Kg 
(at  launch) 


COORDINATE  SYSTEMS 


CJt 


Y*=  Yg  (at  launch) 


FIGURE  1 
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Let  R  be  the  vector  from  the  origin  of  the  system  to  the 

mass  centre  of  the  shell,  and  let  R  have  components  (X  ,  Y  ,  Z  )  in(X  Y  Z  ). 

^  'o  oo'  'ooo' 


d  R 

■g-p  and  has  components 


Then  the  velocity  of  the  mass  centre  is  v 
d  X  d  Y  d  Z 

in  (X^  Y^  Z^  ).  At  the  mass  centre  c  ,  we  construct 


d  t 


d  t 


d  t 


(X^  Y^  Z^  )  which  moves  with  c  but  remains  always  parallel  to  X^  Y^  Z^  . 
This  is  illustrated  in  Fig.  1. 

We  will  use  the  subscripts  1,  2,  3  for  the  components  of  vectors  in 
(X^  Y^  Zq  ).  Then,  since  the  mass  m  of  the  shell  is  fixed,  Eq.  (1)  may  be 
written 


m 


m 


m 


d  X 
_ o 

d  t 

d  Y 
_ o 

d  t 

d  Z 
_ 0^ 

d  t 


= 


-  G2  +  F2 


(3) 


where 


d  X 
_ c 

d  t 


=  X 


d  Y 
_ o 

d  t 
d  Z 

o 
d  t 


=  Y 


=  Z 


(4) 


Now,  principal  axes  and  moments  of  inertia  for  the  shell  may  be  computed. 
The  c-xyz  coordinate  system  may  be  constructed  with  the  y- axis  along  the 
principal  longitudinal  axis  of  inertia  (plai),  and  the  x  and  z  axes  along  the  other 
two  (transverse)  principal  axes  as  shown  in  Fig.  1.  The  orientation  of  c-xyz  is 


6. 


defined  by  the  three  rotations  4>  ,  ^  ,  ip  taken  in  the  indicated  order.  The 

angular  velocity  vector  of  c-xyz  will  be  denoted  by  .  Let  I^,  1^, 
denote  the  principal  moments  of  inertia  along  the  x,  y,  and  z  axes.  It  will  also 

/\  y\  ^ 

be  convenient  to  use  i,  j,  k  as  unit  vectors  in  the  x,  y,  z  directions  respectively. 
We  will  use  the  subscripts  x,  y,  z  for  components  of  vectors  in  c-xyz  .  The 
angular  momentum  of  the  shell  (about  c  ),  can  now  be  written  as: 


/s.  A*  ^ 

ilo)  +  jIoj  +kla) 
XX  y  y  z  z 


(5) 


From  Fig.  1  it  is  apparent  that  the  components  of  angular  velocity  are  given  by: 


OJ 


d  4*  d  6 

cos  6  sin  ip  +  ^  ' 


OJ 


d  4>  a  ,  d  ip 

-T—  Sin  o  +  ^  , 

d  t  d  t 


0J_ 


d  d  6  . 

-r—  cos  9  cos  ip  +  -r—  sm  ip  , 
d  t  d  t  ^ 


which  can  be  inverted  to  vield: 


d  4> 

d  t 


sec  9 


w  cos  ip  -  ui  sin  ip 

Z  X  J 


d  9 
d  t 


CO  cos  ip  +  (jj  sin  ip  , 
X  z 


(6) 


d  ip 
d  t 


=  CO  -  tan  9 


'  CO  cos  ip  -  u  sin  (il 
L  z  ^  X 


(We  assume  that  9  is  always  within  the  bounds  I  9  |  < 

d  i” 

Differentiating  Eq.  (5),  and  recalling  that  = 


rr/2  ). 

_  yv 

CO  X  i,  etc.  we  obtain 


o 
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je 


d  t 


i  f  (I  oo)  +  (I  -I)coco] 
Ldt  X  x'  'z  y  zyj 

+  j  [  -  (I  00)+(I  -I  )C0  CO  1 

•’Ldt  yy  X  z'xzj 

+  k  [  (Iw)+(I  -I)cocol 

Ldt  'zz'  'y  x'yxj 


We  now  assume  that  the  shell  is  rotationally  symmetric  so  that  we  can 

set  1  =1  =  I,  1  =  r.  Then  Eq.  (2)  can  be  written  in  the  scalar  form: 

X  z  y  i  \  ^ 


d  CO 
_ x 

d  t 


(1  -  r)  CO  CO  +  N 
y  z  X 


d  CO 

r  =  N 

d  t  y 

d  CO 


(7) 


I  ^  ^  (I-Duo  + 


N 


When  we  consider  the  aerodynamic  forces  and  moments  it  will  be 
convenient  to  have  the  velocities  relative  to  the  undisturbed  air  in  the  neighborhood 
of  the  shell  rather  than  the  fixed  X  Y  Z  coordinate  system.  Therefore 
in  the  neighborhood  of  the  shell  the  undisturbed  air  flow  would  have  had  a  linear 
velocity  W  .  The  velocity  of  the  shell  relative  to  air  can  then  be  defined  by; 


so  that 


V  =  v  -  W 


X  = 

V, 

+ 

W 

c 

1 

$ 

Y  = 

+ 

w 

o 

2 

Z  = 

Vo 

+ 

w 

o 

3 

(8) 


8. 

where  are  known  functions  of  time  and  position.  We  can  now 

substitute  Eq.  (8)  into  Eq.  (3)  and  (4)  to  obtain 


d  X 
_ o 

d  t 


Vf  +  Wi 


d  t 


^2  +  ^2 


d  Z 
_ o 

d  t 


V3  + W3 


(9) 


and 

d  V 

m  — -3— —  =  “mW,+F,+G, 
d  t  ill 

d  V2 

"d^  =  -mW2  +  F2  +  G2  (10) 

d  V 

.  d  W 

where  we  have  set  =  -g-y-  >  etc. 

The  wind  velocities  will  be  specified  in  a  coordinate  system  fixed  in 

the  earth.  As  mentioned  previously,  they  are  specified  by  a  direction 

(azimuth),  plus  vertical  and  horizontal  components  P.,,  and  Q,,,  .  Let  us  use 

w  w 

the  transformation  of  a  vector  from  earth  to  fixed  coordinates  as  indicated  in  the 


appendix  (A- 23,  A-24). 

Then,  the  wind  velocities  become: 


9. 


Wi  =  (Ag  (A'^  cos  fi  t  +  D’j^  sin  n  t )  +  ('A'^^  sin  t+  cos  n  t) 


+  S  °'l  > 


+  (Gg  (A'^  cos  t  +  D'^  sin  n  t)  +  H^,  (“A'j^  sin  n  t  +  cos  n  t ) 


^2  ^  ^  ^'l  t )  +  Bg  (-B'j^  sin  n  t  +  E'j^  cos  n  t  ) 


^E  L  ^ 


+  Qw  ^^E  ^  ^'l  ^  (-B'j^  sin  n  t  +  E'^  cos  fit) 


■^Je»'l> 


^3  ^  ^W  ^■^E  ^  ^  L  ^  )  +  ®E  L  ^  cos  n  t  ) 


+  Q^  (Gg  (C'j^  cos  n  t  +  F’j^  sin  n  t  )  +  Hg  (-C^  sin  12 1  +  cos  S2  t  ) 


^E  L  ^ 


The  quantities  A^,  . . .  J^,  A'^  . . .  J'^  arise  due  to  the  transformation 
as  discussed  in  the  appendix,  and  are  defined  there. 

Eqs.  (6),  (7),  (9)  and  (10)  constitute  a  system  of  twelve  first  order 
differential  equations  in  the  twelve  dependent  variables  X  ,  Y  ,  Z  ,  V,,  V„,  V„, 
>  9  I  ^  ,  CAjy,  ,  and  the  independent  variable  t  (time  ).  It  remains 

now  to  express  the  other  terms  in  this  set  of  equations  in  terms  of  these  twelve 
generalized  coordinates  and  the  time. 
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We  will  need  expressions  for  the  direction  cosines  of  the  X  Y  Z 

o  o  o 

and  c-xyz  axes  with  respect  to  each  other.  From  Fig.  1  to  within  first  order 
terms  (in  4> )  we  have: 

=  cos  ip  -  ^  sin  6  sin  ip  , 

^x2  ~  ^  ^  ^  ^  ^  ’ 

a  „  =  -cos  0  sin  ip  , 
xo 


^2 

^3 


-4)COS  0  , 
cos  0  , 
sin  0  , 


(12) 


) 


a  =  sin  ip  +  sin  0  cos  ip  , 

Z  li 


a^2  "  ‘t’  ^  ~  ®  ^  > 


a^2  =  cos  0  cos  ip  . 


Then  if  u  is  a  vector  we  have 


^x  =  ^xl  ^1  +  ^x2  ^2  +  \3  ^3  ■ 


u,  =  a,u  +a,u+a,u  ,  etc. 
1  xi  X  yi  y  zi  z 


! 
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a)  Those  independent  of  angular  velocity  (Drag  Force,  Lift  Force 
or  Force  due  to  cross  velocity,  Restoring  Moment  or  Moment  due 
to  cross  velocity), 

b)  Those  depending  on  transverse  angular  velocity  of  the  shell 
(Force  due  to  cross  spin.  Damping  Moment  or  Moment  due  to  cross 
spin), 

c)  Those  depending  on  the  axial  spin  of  the  shell  (Roll  Moment  or  Moment 
due  to  spin,  Magnus  force  and  moment  due  to  cross  velocity  and  cross 
spin). 

Although  these  forces  and  moments  interact  with  each  other,  they  are  traditionally 
treated  separately. 

It  is  also  customary  to  describe  the  aerodynamic  forces  and  moments  by 
means  of  dimensionless  aerodynamic  coefficients.  For  this  purpose  we  introduce 
the  average  mass  density  of  the  air  surrounding  the  shell  and  denote  this  density 
by  p  ;  we  also  need  a  length  characteristic  of  the  shell  and  therefore  introduce 
the  diameter  d  of  the  shell.  Also  we  use  the  linear  and  angular  velocities  of 
the  shell  where  appropriate.  Since  we  wish  to  base  our  computation  on  motion 
of  the  shell  in  undisturbed  air  we  must  account  for  winds  by  using  V  as  the 
velocity  of  the  mass  centre  of  the  shell  with  respect  to  the  air. 

Drag  F'orce  and  Lift  Force.  The  aerodynamic  force  which  acts  on  the  shell  when 
there  is  no  angular  velocity  (w  =  0  )  is  resolved  into  two  components;  the  drag 
force  D  which  is  taken  parallel  to  V  ,  and  the  lift  force  L  perpendicular  to 
V  .  We  define  the  drag  coefficient  Kj^  and  the  lift  coefficient  by 


13. 


From  physical  considerations  we  can  argue  that  there  must  be  at  least  three 
principal  directions  such  that  if  V  is  parallel  to  one  of  them  the  lift  force  is 
zero.  The  principal  direction  most  nearly  parallel  to  the  plai  of  the  shell  is 
called  the  aerodynamic  axis  and  the  angle  between  the  aerodynamic  axis  and  V 
is  called  the  angle  of  attack  (denoted  by  a  ).  Let  ^  be  a  unit  vector  along  the 
aerodynamic  axis.  We  assume  that  a  is  a  first  order  quantity  and  use  the 
angles  and  to  orient  e  in  the  c-xyz  coordinate  system  as  shown 
In  Fig.  2.  Then  to  first  order  terms 


FIGURE  2 
AERODYNAMIC  AXIS 
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Also 


2  .2  .2 

a  =  £  +  £ 

^  X  ^  z 


We  assume  that  the  lift  force  vector  lies  in  the  plane  defined  by  e  and 
V  ,  and  since  we  have  specified  L  perpendicular  to  V  we  could  write 


^  -  K,  pd  Vx(exV) 


=  p  d^  [  ^  -  V  .  e  V  ] 


where 


to  first  order  terms.  We  could  also  write  V  -  g-  =  V  cos  a  =  V  to  first  order 
terms  so  that 

Lf  =  p  d^  V  [  V  (-  4.  cos  e  +  sin  cos  J  , 

^2  "  P  V  [  (V  cos  e  -  V2)  -  V  cos  ip  sin  ^  )  sin  g]  (15) 
L3  =  p  d2  V  [  (V  sin  0  -  V3 )  +  V  cos  ip  +  ^^sin>p)  cos  g] 

The  drag  force  is  written  as 


D  =  -Kj,  p  V  V  , 

SO  that 

Df  =  -Kj^  p  d^  V  , 

D2  =  -Kj3  p  d^  V  V2  . 
D3  =  -Kj^  p  d2  V  V3  . 


(16) 


o 


For  small  angles  of  attack  both  and  can  be  taken  as  independent 
of  a  since  both  must  be  (approximately)  even  functions  of  a  and  therefore  the 
dependence  will  be  second  order.  The  other  significant  variables,  chief  among 
which  is  the  velocity,  are  assumed  to  enter  only  by  way  of  the  Mach  number 
related  to  V  and  we  therefore  allow  and  to  be  functions  of  Mach 

number.  If  we  further  assume  that  both  the  sonic  velocity  and  the  mass  density 
p  depend  only  on  altitude,  the  drag  and  lift  forces  can  be  determined  from  the 
generalized  coordinates  and  the  time. 

Restoring  Moment.  The  aerodynamic  moment  which  acts  on  the  shell  when  there 
is  no  angular  velocity  has  a  component  perpendicular  to  the  aerodynamic  axis. 
Since  this  moment  is  due  to  the  same  forces  as  the  Lift  and  Drag  force,  it  follows 
that  the  restoring  moment  must  be  perpendicular  to  V  also.  We  use  P  for  the 
restoring  moment  and  write 

P  =  Kp  p  d^  , 

where  Kp  is  the  restoring  moment  coefficient,  positive  if  P  •  'e  x  V  ^  0,  It 

is  usually  more  convenient  to  define  a  centre  of  lift  by  the  vector  from  the 

•  r 


mass  centre. 


such  that 


j^p  =  -jep^ 


P  =  £p  X  L  , 


so  that  ip  has  the  same  sign  as  Kp  ,  It  follows  then  that 


P  =  K,  P  i  X  [  V  X  ( 


e  X  V 


=  K.  p  d^  ip  (e  .  V  )  (e  X  V  )  . 


16. 


and  to  first  order  terms 

r»  _ 

^  ■^P  ^  L  ^^1  ^  ^2  ^  ^  ~  ^^2  ®  ~  ^3  ®  ^ 

-  ^  ^3  ®  )] 

^y  "  °  ’  (17) 

^z  "  '^L  P  -^P  '^  [  <’^1  +  "^2  ^  +  (V^  sin  0  -  Vg  cos  0  )  sin  ^ 

+  (V2  cos  6  +  Vg  sin  0  )J 

It  is  presumed  that  ^p  is  a  known  function  of  time  and  the  Mach  number  associated 
with  V  .  Dependence  on  all  other  parameters,  including  angle  of  attack, is 
ignored.  Generally,  the  location  of  the  centre  of  lift  is  taken  as  a  function  of 
Mach  number  only,  while  the  location  of  the  mass  centre  is  a  function  of  time  only. 
Force  and  Moment  Due  to  Cross  Spin.  If  the  shell  has  an  angular  velocity  w 
the  resultant  aerodynamic  force  and  moment  will,  in  general,  be  different  from  the 
case  when  angular  velocity  is  absent.  For  if  p  is  a  vector  from  the  mass  centre 
to  some  point  on  the  surface  of  the  shell,  the  velocity  of  this  point  will  be  V  +  w  x 
with  respect  to  the  air,  and  among  other  things,  the  aerodynamic  force  exerted  on 
a  small  part  of  the  surface  depends  on  the  velocity  of  that  surface  with  respect  to 
the  nearby  undisturbed  air.  The  component  of  the  angular  velocity  which  is 
perpendicular  to  the  aerodynamic  axis,  ,  is  called  the  cross  spin.  Clearly 

Wg-ex(coxe)=u)-(^-u))'e. 

The  additional  aerodynamic  force  due  to  cross  spin  is  denoted  by  S  and  is  assumed 
to  be  perpendicular  to  both  ^  and  w  ; 
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Q 

S  =  K  pd  Vex” 
s  ^  s 

=  K  p  V  ^  X  oj 


where  is  the  cross  spin  force  coefficient.  To  within  first  order  terms, 

3 

Si  =  -  Kg  p  d  V  ^  )  sin  !/»  -  )  cos 

3  r 

^2  ""  ^s  P  ^  ^  [  ^“x  “y  ^  ip+  (<^2  ~  ^x  “y  ^  ® 

^3  ^  “^s  P  ^  [  ^“x  ■'■  “y  ^  4>  +  “y  )  ® 

A  similar  argument  can  be  advanced  for  the  existence  of  a  moment  due  to 
cross  spin,  H  ,  which  must  be  parallel  to  the  cross  spin.  We  therefore  write 

H  =  -K„  p  d^  V  w 

rl  S 

'  "Kh  P^i^Ve  X  (w  X  e) 
so  that,  to  first  order  terms, 

Hx  =  -K^  p  d^  V  [u,^  +  u;y] 

“y  "  °  ’  (19) 

“z=  -^hP^^v 


where  Kj^  ,  the  cross  spin  moment  coefficient  is  also  called  the  damping  moment 

coefficient.  This  name  follows  from  the  fact  that  for  p'  0,  the  moment  H 

rl 

is  proportional  to  the  cross  spin  and  tends  to  damp  it  out.  ^ 

The  coefficients  Kg  and  are  usually  taken  to  equal  zero  since  their 
effect  is  generally  small  compared  to  the  effects  of  other  causes.  In  any  event. 
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we  presume  that  and  are  at  most  known  functions  of  the  Mach 

number  associated  with  V  .  Sometimes  it  is  convenient  to  define  a  cross  spin 

centre  of  pressure  by  the  vector  I  from  the  mass  centre,  £  =  ^  where 

o  s  s 

Kpj  =  Kg  £g  .  Note  that  £g  can  now  be  a  function  of  Mach  number  and  time. 

Also,  if  the  cross  spin  centte  of  pressure  is  behind  the  mass  centre  (£  >  0  ) 

s 

the  damping  moment  does  indeed  oppose  the  cross  spin. 

Roll  Moment.  To  this  point  we  have  accounted  for  the  aerodynamic  force  and 
moment  due  to  the  linear  and  angular  velocity  of  the  shell  in  still  air,  except  for 
the  component  of  moment  along  the  aerodynamic  axis  due  to  linear  velocity,  and 
the  force  and  moment  due  to  the  component  of  angular  velocity  along  the  aero¬ 
dynamic  axis.  Aerodynamically,  we  assume  that  the  shell  is  a  surface  of 
revolution  about  the  aerodynamic  axis.  Due  to  skin  friction,  the  component 
of  angular  velocity  along  the  aerodynamic  axis  causes  a  friction  moment  about 
this  axis.  If  R  denotes  this  moment,  then  we  can  write 

R=  -Kp  p  d^V(  «.e)e  , 

or  to  first  order  terms, 

%  =  -Kp  p ,  (20) 

= -Kp  p  dS  • 

We  suppose  that  Kp  ,  the  roll  friction  coefficient,  is  at  most  a  function  of  the 
Mach  numbers  associated  with  V  and  ^  w  . 
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Magnus  Force  and  Moment.  We  have  finally  to  consider  the  Magnus  force  and 
moment.  These  are  significant  only  when  there  is  appreciable  angular  velocity 
along  the  aerodynamic  axis.  Then  there  will  be  a  force  perpendicular  to  both 

—  As 

V  and  e  (if  they  are  not  parallel  )  which  can  be  represented  by 

P  o)^  X  V 
MM  e 

where  =  (  w  •  e  )  "e  is  the  component  of  angular  velocity  along  the  aero¬ 
dynamic  axis,  and  Kj^  is  the  Magnus  force  coefficient  due  to  cross  velocity, 

assumed  to  be  at  most  a  function  of  the  Mach  number  associated  with  —  to  and 

2  e 

V  .  To  first  order  terms 

^M  =  ^M  P  ‘^y  ^  ^  V 

so  that  * 

^M1  "  '^^M  P  “y  [  ^  ^^x  sin  ip)  +V2  sin  0  -  V3  cos  e] 

3  r 

^M2  ■  ^M  °  “y  L  "^1  ®  ^3  ^  ®  '  ^x  ^  ^  ^  ] 

3  r 

^M3  ’  '^M  ^  “y  L^l  ®  V2  (  4>  cos  0  -  sin  cos  p  )] 

Now  because  the  cross  spin  adds  a  cross  velocity  to  poiits  away  from 

the  mass  centre  (in  the  amount  to^  x  p  where  p  is  the  position  vector  from 

the  mass  centre  )  we  might  expect  the  cross  velocity  to  affect  the  Magnus  force. 

We  therefore  suppose  Kj^  to  be  measured  at  zero  cross  spin  and  consider  a 

Magnus  force  due  to  cross  spin,  F./  with  coefficient  .  We  therefore 

ivi  M 


(21) 


write 


20. 


— 


F\,  =  K'.,  pduex(exw) 
M  M  y  '  s 


4  _ 

=  -K\,  p  d  CO  oj 
M  y  s 


(21) 


-  -K'mP"  ‘^y 


[^<“x 


+  CJ  ^  )+k(u  -CO  £  ) 

y  z  z  y  ^x ' 


so  that 


p’mi  = 


,  p  d^  CO  [  (co  +  CO  £  )  cos  0  +  (cj  -  CO  £  )  sin  ibl 
M  y  L  X  y^z'  ^  '  z  y  '’x ' 


F'  _  =  -K'  ^ 
M2  M 


P  ^  Wy  "^^y  ^  ^  ^'^z  '^y  ^x  ^  .]  sin  0  ,  (22) 


^'m3  "  *^'m  P  “y  [^“x  “y  ^z  ^  >P  '  ^  ®  ■ 


We  should  also  expect  a  Magnus  moment  due  to  cross  velocity  and  one  due 
to  cross  spin.  We  define  a  Magnus  cross  velocity  centre  of  pressure  by  the  vector 
J^M  ~  ’^M  ^  the  mass  centre  and  allow  to  be  a  function  of  the  Mach 

number  associated  with  V,  ^  oo^  and  possibly  the  time.  Then  we  write 


-  -K^  P  “y  '  *  > 

=  %  )®] 

SO  that 

“mx  =''MP‘''«M“y  +(V^ +V^<|)  )  cos^  +  (V^  sin  6  -  Vg  cos 0  )  sin ^1 

M^y  =  »  .  (23) 

“mx  -  ''M'’‘‘^V“y 
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Similarly,  we  define  a  Magnus  cross  spin  centre  of  pressure  by 

,*6  and  write 
M  M 


M 


X 

F  ' 

M 

M  P 

'*  ^  M 

0) 

y 

e 

X  CO 

s 

.4  . 

[a 

M  P 

ca 

y 

i  (co  -  CO 
-  Z  y 

y 


0  "T  Cl) 

X  y  'Z 


so  that 


M' 


Mx 


-  p  d  -  w  ,  )  , 


y  '  z  y  ^x 


M’My  =°’ 


M' 


Mz 


M  y  '  X  y  ^z 


(24) 


IV.  SUMMARY 
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The  equations  of  motion  along  with  the  necessary  notation  are  collected 
in  this  chapter  for  easy  reference.  The  independent  variable  is  time  denoted  by  t  . 
The  equations  of  motion  are  in  the  form  of  a  system  of  twelve  first  order  differential 
equations  in  twelve  coordinates: 


X„Y  Z 
o  o  o 


Vi  V3 


<t>.  0.  ^ 


components  of  the  displacement  vector  of  the  mass  centre  of  the 
shell  from  the  launching  site  (measured  in  inertial  coordinates) 
inertial  components  of  the  velocity  vector  of  the  mass  centre 
relative  to  the  air 
angles  orientating  the  shell 

components  of  the  angular  velocity  vector  of  the  shell 


The  twelve  equations  of  motion  are 


d  X 

0 

d  t 

d  Y 

o 
d  t 

d  Z 
_ o 

d  t 


=  Vi+Wi 


V2+W2 


^3  +  ^3 


d 

m  -g-p  =  -mW^  +  F^  +  G^ 
dV2 

m  _  =  -mW2  +  F2  +  G2 

d  V3 

-dT  =  -"^Wg+Fg+Gg 


(9) 


(10) 


m 
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d  <)) 

irr 

=  sec  ®  ^  “2  ^  ”  ‘^x  ^  ] 

d  e 

d  t 

=  cos  4>  +  ^2  ^ 

(6) 

d  jp 
d  t 

=  Wy  -  tan  ®  [  cos  ^  sin  ^  j 

I  4^  =  (r-i)<.y<.^+N^ 


/ 

I 


N 

y 


^  -d-r  = 


(7) 


We  now  define  the  individual  terms  in  the  above  equations: 


t  +  D'^  sin  fit)  +  B^.  (-A’^  sin  n  t  +  D’^  cos  fit) 


+  Ce  G-^) 


+  Qw  L  52  t  +  sin  n  t  )  +  Hg  (-A'^  sin  q  t  +  D'^  cos  n  t  ) 


■'JeG’l> 


^2  ^^E  L  52  t  +  E'j^  sin  n  t  )  +  (  -B' sin  n  t  +  E'^^  cos  n  t  ) 


+  ‘^e^’l> 


(11) 


+  Qw  ^‘^E  ^  t  +  E'^  sin  n  t  )  +  Hg  ( -B'^  sin  n  t  +  E' cos  n  t  ) 
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d  4) 
d  t 

■  =  sec  6  [  ‘*^2:  ^  ~  ‘^x  ^  J 

d  0 
d  t 

=  cos  ip  sin  ip 

(6) 

d  ill 
d  t 

=  ojy  -  tan  ®  cos  ^  sin  ^  j 

d  t 

.  d  oj 


(r  - 1)  oj  oj  +  N 
y  z  X 


=  N 


(7) 


d  t 


(I  -  r)  oj  oj  +  N 

y  X  z 


We  now  define  the  individual  terms  in  the  above  equations: 


''"l  ■  t  +  D’l  sin  R  t  )  +  Bg  (-A’^  sin  R  t  +  D’^  cos  R  t  ) 


+  Qw  <Ge  <A'l  cos  R  t  +  D’g  sin  R  t  )  +  Hg  (-A'g  sin  R  t  +  D’g  cos  R  t  ) 


-MeG’l) 


^2  ^^E  R  t  +  E’g  sin  R  t  )  +  Bg  (-B'g  sin  R  t  +  E’^  cos  R  t  ) 


+  CeA’l) 


(11) 


+  Qw  ^^E  +  E’^  sin  R  t  )  +  Hg  (-B’^  sin  Rt  +  E’^  cos  R  t ) 


+  Je  »’l> 
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^3  "V  +  F’l  sin  n  t )  +Bg(-C'j^  sin  J2 1  +  F' 


+  S  J’l  > 


■•■^W  L  ^°s  +  F'l  sin  n  t)  +  H  (-C’l  sin  J2  t  +  F’ 


+  JeJ’l> 


Pq  2 

Fl  =  -m  (  —  )  (  Gg  (A'^  cos  n  t  +  sin  fi  t  ) 

+  Hg  (-A'j^  sin  n  t  +  D*g  cos  fi  t  )  +  Jg  G’  ) 


o  .2 


^2  ■  ^o  ^  ^  ^  c°s  fi  t  +  E'g  sin  fi  t  ) 


Pe 


+  Hg  (-B'g  sin  n  t  +  E'g  cos  fi  t  )  +  Jg  H'g  ) 

Po  2 

^3  =  So  <7“)  (  C3g  (C'g  cos  n  t  +  F'g  sin  n  t  ) 
+  ^^(-C'g  sin  n  t  4-  F'g  cos  52  t  )  +  Jg  J'g  ) 


Ag  -  cos  6  cos  <t)g  cos  iji  +  sin  0p  sin  ip^ 


Bg  =  sin  0  cos  <()  cos  ip  -  cos  0  sin  ip 


Cg  =  -sin  <1^  cos  i/jg 


Gg  =  cos  0g  sin  <t)g 


Hg  =  sin  0g  sin  <t)g 


cos  n  t ) 

cos  n  t ) 

(11) 


(13) 


(A -25) 


Jg  =cos<t)g 


c 
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Pe^  =  V  +  + 


‘t>E  =  Zg  /  +  Yg-  +  Zj. 


i  -  '  I-  2  .,.  2  _  2 


Og  =  sin 


Yg/  Jxj+Y.^ 


E  ^ 


^E  ^o  ^  ^L  ^  t )  +  Y^  (Dg  cos  n  t  +  Eg  sin  fi  t ) 

+  (Zq  +  p  )  (Gg  cos  n  t  +  Hg  sin  fi  t  ) 

'^E  ®i"J2t  +  Bg  cos  nt)+Y^(-Dg  sin  fi  t  +  Eg  cos  i2  t ) 

+  (Zq  +  p  )  (“Gg  sin  n  t  +  Hg  cos  n  t ) 


Ag  =  cos  Sg  COS  4)g  cos  ^g  +  Sin  0g  sin  ip^ 
Bg  =sin  Og  cos  cos  ip^  -  cos  0g  sin  ip^ 
Cg  =  -  sin  4>g  cos  i^g 

Dg  =  cos  0g  cos  4)g  sin  ^g  -  sin  0g  cos  ip^ 
Eg  =  sin  0g  cos  4)g  sin  ^g  +  cos  0g  cos  ip^ 
F  =  -  sin  4)  sin  ip, 

•L»  Li  L 


Gj  -  cos  0  sin  cf) 

•Lj  Li  E 


H-  =  sin  0.  sin  4) 


L  X 


Jg  =COS4>g 


(A-19) 


(A -17) 


(A- 18) 
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A’^=  cos  cos  <))^  cos  +  sin  ip^  sin  0^ 

B'j^  =  sin  ip^  cos  cos  0^^  -  cos  jp^ 


C’j^=  sin  <l)j^  cos  0j^ 


=  cos  ip^  cos  <|)^  sin  0^^  -  sin  cos  0^^ 
E'j^=  sin  cos  sin  0^^  +  cos  ipj^  cos  0^^ 


(A- 


F’^=  sin  sin  0^ 


G'.  -  -cos  ip  sin  <(> 


H’  =  -sin  >p  sin  <(> 


J'l  =  cos 


Gi  =  D,  +  S  +  +P'm1 


=  -K^  p  V  V.^  , 


D2  =  -Kj^  P  d  V  V. 


(16) 


D3  =  -Kj3  p  d  V  V3 


K  "  ^  V  (  V  (-<t>  cos  0  -I-  sin  }p  -  cos  ip  )  -  Y ^  ) 

2 

1^2  =  p  d  V  (  (V  cos  0  -  V2  )  -  V  cos  ^  sin  )  si  n  0  )  (15) 

2 


13) 


4 
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28. 


R 

=  K„  p  d^  V  00  ^ 

X 

F  y  z 

R 

y 

=  "K,:.  p  d^  V  00 

F  y 

(20) 

R 

=  -K„  p  d^  V  00  ^ 

z 

F  ^  y  ^x 

^Mx  «M“y<  V  +  (V^  +  4>  )  cos  ip  +  (V^  sin  ©  -  cos  0)  sin;/)) 

^My  =  ^  (23) 

^Mz  "  ^M‘^y  ^  ^  ‘  ^^3  ‘^osB)  cos  ify 

+k;.  .  p  d^  i\.  00  (w  -  CO  e  )  , 

MX  M  My  z  y  x 


MW  \W  i  w  u 

y  X  y 


(24) 
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Notation 


+  V 


2 

3 


W^,  W^,  Wg  -  components  of  the  wind  velocity  in  the  neighborhood  of  the  shell, 
relative  to  the  inertial  reference  system 


m 

I 

I’ 


P 
d 

Df,  D^,  Dg 

^1’  S’  S 

Sf,  S2,  Sg 


Ml’ 


M2’  "  M3 


^  Ml’  ^  M2’  ^'m3 

P  .  P 
X  y  z 

»x’  S’  Hz 
S’  S’  Rz 

M,,  ,  M-.  ,  M 
Mx  My’  Mz 

,  M\,  ,  M’ 

Mx  My’  Mz 


K 


D 


K, 


K. 


K 


M 


mass  of  the  shell  (constant) 

principal  transverse  moment  of  inertia  of  the  shell  (constant) 

principal  longitudinal  moment  of  inertia  of  the  shell  (constant) 

aerodynamic  malalignment  angles  (constant) 

density  of  the  air  (function  of  position),  i.  e.  Oj-,,  p  <}) 

shell  diameter 

aerodynamic  drag  force 

aerodynamic  lift  force 

aerodynamic  cross  spin  force 

Magnus  force  due  to  cross  velocity 

Magnus  force  due  to  cross  spin 

aerodynamic  restoring  moment  (due  to  cross  velocity) 
aerodynamic  damping  moment  (due  to  cross  spin) 
aerodynamic  roll  moment 
Magnus  moment  due  to  cross  velocity 
Magnus  moment  due  to  cross  spin 

aerodynamic  drag  coefficient  (function  of  Mach  number) 
aerodynamic  lift  coefficient  (function  of  Mach  number) 
aerodynamic  cross  spin  force  coefficient  (function  of  Mach  number) 
Magnus  force  coefficient  due  to  cross  velocity  (function  of  Mach 


number) 
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o 


^w’ 


Pe 

©L’  ’  Pl 


n 


aerodynamic  damping  moment  coefficient  (function  of 
Mach  number) 

aerodynamic  roll  fricticn  coefficient  (function  of  Mach  number) 
distance  between  mass  centre  and  centre  of  lift  (function  of 
Mach  number) 

distance  between  mass  centre  and  Magnus  cross  spin  centre 
of  pressure  (function  of  Mach  number) 
horizontal  and  vertical  components  of  wind  velocity 
direction  of  wind  velocity  as  given  in  coordinates  fixed  in 
the  earth 

coordinates  of  a  point  measured  in  earth  coordinates 
coordinates  of  the  shell  at  time  of  launch  (fixed  with  respect 
to  the  earth) 

angular  velocity  of  the  earth  =  2  it  radians/day 
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K 


H 


^w’  ‘^w 

®E’*E’  Pe 
®L’  ‘*’l  ’  Pl 


J2 


aerodynamic  damping  moment  coefficient  (function  of 
Mach  number) 

aerodynamic  roll  frictioi  coefficient  (function  of  Mach  number) 
distance  between  mass  centre  and  centre  of  lift  (function  of 
Mach  number) 

distance  between  mass  centre  and  Magnus  cross  spin  centre 
of  pressure  (function  of  Mach  number) 
horizontal  and  vertical  components  of  wind  velocity 
direction  of  wind  velocity  as  given  in  coordinates  fixed  in 
the  earth 

coordinates  of  a  point  measured  in  earth  coordinates 
coordinates  of  the  shell  at  time  of  launch  (fixed  with  respect 
to  the  earth) 

angular  velocity  of  the  earth  =  2  it  radians/day 
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V.  APPENDIX 

COORDINATE  SYSTEMS 

i  We  will  first  derive  some  general  transformations,  and  then  apply  them 

to  specific  applications.  Let  us  consider  4  such  transformations. 


ii  General  Transformations 


ii-1  Rotation  of  Coordinate.  Axes 


Let  us  successively  make  rotations  about  the  X,  Y,  and  Z  axes  res¬ 
pectively;  through  the  angles  9,  <t>,  and  tp  .  Then  we  will  obtain  the  following 
relations  between  unit  vectors: 


A'  A  A  jA 

i'=  iA+jB+kC 


i  =  i'  A'  +  j'  B'  +  k'  C 


A  ^  A 

j'=iD+jE+kF 


A 

j  =  i'  D'  +j'  E'  +  k'  F' 


A-1 


^  A  A  A 

k'=  iG-l-  j  H  -l-kj 


✓V  ^  A 

k=  i'  G'  +  j'  H'  +  k'  J' 


A  A  ^  ^  Av. 

Here  i,  j,  k,  and  i',  j',  k',  are  unit  vectors  along  the  X,  Y,  Z,  and  X',  Y',  Z' 


axes  respectively. 


Also: 

A  =  cos  0  cos  4>  cos  ip  +  sin  6  sin  ip 
B  =  sin  9  cos  <t)  cos  ip  -  cos  9  sin  ip 
C  =  -sin  4)  cos  ip 

D  =  cos  9  cos  4>  sin  ip  -  sin  9  cos  ip 
E  =  sin  9  cos  4>  sin  ip  +  cos  9  cos  ip 
F  =  -sin<))  sin  ip 
G  =  cos  9  sin  <!> 

H  =  sin  9  sin  <)) 

J  =  cos  4) 


Figure  A-1 

A’  =  cos  ip  cos  4>  cos  9  +  sin  ip  sin  9 
B’  =  sin  ip  cos  4>  cos  9  -  cos  ip  sin  9 
C  =  sin  4’  cos  9 

D'  =  cos  ip  cos  sin  9  -  sin  ip  cos  9 

E'  =  sin  ip  cos  4>  sin  9  +  cos  ip  cos  9 
F'  =  sin  6  sin  9 
G'  =  -cos  ip  sin  4* 

H'  =  -sin  ip  sin  4> 

J'  =  cos  4> 


Then  any  point  may  be  located  by  the  vector: 
7  =  iX  +  j'Y  +  kZ  =  i'X’+r  Y'  +  k’  Z’ 


34, 


where; 


X  =  AX'  +DY'  +  GZ' 

X’  =  A'X+  D’Y+  G'Z 

Y  =  BX'  +  EY'  +  HZ’ 

Y'  =  B'X  +  E'Y  +  H’Z 

A- 3 

Z  =  CX'  +  FY'+ JZ' 

Z’  =  C'X+  F'Y  + J'Z 

ii-2  Translation  of  Coordinate  Axes: 


Consider  X  Y  Z  to  be  a  rectangular  system  of  coordinates  fixed  in  the  earth, 
with  the  Z  axis  along  the  axis  of  rotation  of  the  earth.  Let  X*,  ^  be  an 

inertial  system  such  that  Z+  will  coincide  with  Z  .  The  X,  Y,  Z  system  then  rotates 

A 

with  a  constant  angular  velocity  k  S2  with  respect  to  the  X*  Y^  Z^  system. 


n=  2*Tr 

day 


1  ^  cos  t  +  ^  sin  f2  t; 
j  =  -i*  sin  f2  t  4  cos  Q  t 
k  =  k 


A  ^ 

=  i  cos  n  t  -  j  sin  n  t 

A  ^  ^ 

4  =  i  sin  n  t  +  j  cos  n  t 


k 


A, 

k 


Then,  any  vector  may  be  represented  by: 

r-iX+jY4-kZ  =  i  X  +  j  Y  +  k  Z 

*  *  •'  *  *  +  * 

where’ 

X  =  X^  cos  r2  t  +  Y^  sin  n  t;  X^  =  X  cos  f2 1  -  Y  sin  f2 1 

^  ~  *^  ^  '*'  ^*  ^  ^  f2  t  +  Y  cos  n  t 

ii-4  Spherical  Coordinates; 


A-5 


A-6 


It  will  be  convenient  to  use  spherical  coordinates  in  addition  to  the  rectangular 


coordinates. 


I  p  sin  9  sin  6 


p  sin  <)>  cos  6 


We  have: 


F igure  A- 4 


X  =  p  sin  4)  cos  0 
Y  =  p  sin  4)  sin  6 
Z  =  p  cos  4> 


2_  ^2  ^  +  2? 


4  =  cos"  Z/7x  +  Y  +  Z 


e=  sin"^  Y//x^  +  Y^ 


Location  of  Points;  and  Direction  of  Vectors  in  Space 


Pertinent  quantities,  such  as  air  densities,  wind  velocities,  and  gravity  forces 
are  available  in  a  reference  system  fixed  with  respect  to  the  earth.  Any  point  in  this 
system  may  be  located  by  specifying  two  angles  0^,  4^  ,  and  a  length  p^  .  These 
quantities  are  the  longitude,  latitude,  and  distance  from  the  center  of  the  earth 
respectively. 

We  may  also  determine  a  direction  at  the  point  (0_  ,  4„  ,  p  ). 

E  E  ^  * 
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"ES' 


ES’ 


Consider  an  additional  rotation  tp^  as  performed  in  section  II- 1.  Then  if 

chosen  so  that  the  direction  (or  vector)  lies  in  the  plane  Y'  0  =  0-  the  vector 

ES 

may  be  represented  by  2  components;  one  along  the  XVc  axis,  and  one  alone  the  Z’ 

tib  &  gg 

axis.  The  angle  4,^  +  180°  is  referred  to  as  the  azimuth  of  the  vector;  and  the  angle  a 
between  the  vector  and  the  Xgg  axis  is  called  its  elevation. 

The  angle  0^  varies  from  0  to  360°.  As  0^  goes  through  the  variations; 

0-^90  —^180  — ^  270  -*-360°  longitude  goes  through  the  variations; 

0-^  90°E  ^180°E  ^  90°W  —0°  . 

The  angle  4)g  varies  from  0  to  180°.  As  (j)^  goes  through  the  variations; 

^90  >-  135  *-180°  latitude  goes  through  the  variations; 

90°N-^45°N— ^0°-^45°S-^0°S  . 
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The  angle  tp^  varies  from  0°  to  360°  .  As  goes  through  the  variation; 
0 — *-90 — »-180 — ►-270*^ — *-360°  ,  azimuth  goes  through 

— *-270^^ — >“360^^ — !fe-90° — *-180°  ,  and  the  more  familiar  directions  go 
through  S  — *“W  -^N  -*-E  — 

iv  Transformation  of  the  Coordinates  of  a  Point  From  a  System  Fixed  in  the  Earth 
t°  the  X’^g  ,  Y’^g  ,  Z’^g  (or  X^Y^Z^)  System 

The  transformation  of  coordinates  from  the  system  fixed  in  the  earth  to  the 
inertial  system  may  be  carried  out  using  4  successive  transformations.  (We  will 
assume  that  the  location  of  a  point  is  specified  by  <J),  0  ,  p  ). 

1.  From  the  (p^  ,  0^  )  system  to  the  X^,  ,  Y^,  ,  system  (both  fixed 

in  the  earth),  using  ii-4. 

Xg  =  pg  sin  4>£  cos  0^ 

Yg  =  pg  sin  sin  0^  A-8 

Ze  =  Pe  °°®  ‘*’e 

2.  From  the  Xg.Yg  Zg  system  fixed  in  the  earth  to  the  X*  Y^ Z  ^  inertial 
system  using  ii-3 

=  Xg  cos  n  t  -  Yg  sin  n  t 

=  Xg  sin  n  t  +  Yg  cos  J2 1  A-9 

Z..  =  Z„ 

3.  From  the  X^Y^Z^  inertial  system  to  the  X'^  ^'o  ^'o  system 

using  ii-1 

+  E',  Y,  +H',  Z 
o  E  p  L  r  L  * 


A-10 


3 

Here  the  angles  appearing  in  A'^ - are  the  angles 

giving  initial  position  and  orientation  of  the  missile  at  launch. 


4.  From  the  X'^,  Y’  ,  Z’^  inertial  system  to  the  X '  ,  Y'  ,  Z'  (or 

Ou  oS  oS 

Y^,  Z^)  inertial  system  using  ii-2. 


Then  the  transformation  becomes; 
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x/ pg  -  sin  (A’j^  (cos  Gg  cos  J2  t  -  sin  6^  sin  J2  t)  +  D'j^  (cos  0^  sin  S 

+  G'j^  cos  <j)g 

y/pg  =  sin  <|)g  (B'^  (cos  0^  cos  J2  t  -  sin  0^  sin  fi  t)  +  E'l(cos  0g  sin  fi  t 

+  H'.  cos  <i) 

L 

z/  pg  =  sin  <l>g(C'g  (cos  0g  cos  J2  t  -  sin  0g  sin  fit)  +  F'^  (cos  0g  sin  fi  t 

+  J'g  cos  (Jjg  -  1 

where: 


-  cos  !p  cos  (h  cos  0 
i-<  J-i  1  1 

+  sin  ip^ 

.  sin  0 

L  L 

«’l  = 

'  sin  0g  cos  41^  cos  0g  - 

cos  ip 

L 

sin  0 

J^ 

11 

b 

sin  4»  cos  0 

II 

Q 

cos  0g  cos  4)g  sin,0g  - 

sin  0g 

cos  0 

J-i 

E» 

^  L‘ 

sin  ipj  cos  4)  sin  0,  + 

cos 

cos  0 

J^ 

II 

sin  4)  sin  0 

J-i  L 

o 

t- ' 

II 

-cos  0  sin  (Jj, 

J-*  J_i 

II 

X 

-sin  ip  sin  4>j 

J^  i_« 

J’l  = 

cos  «► 

J^ 

Similarly,  by  reversing  the  order  of  the  transformations,  we  may 
coordinates  pg,  0g,  4,^  in  terms  of  X^Y^and  . 

5.  Applying  ii-2 


X'  =  X'  „  =  X 
o  oS  o 

Y'o  =  =  Y,, 

^'o  =  ^'oS+  Pe  P] 


!  t  +  sin  0g  cos  fit)) 

+  sin  0g  cos  fit)) 
A-12 

+  sin  0g  cos  fit)) 


A-13 


express  the 


A-14 


6.  Applying  ii-1 


X 

A,  X'  +  D,  Y'  +  G, 

Z' 

* 

L  o  L  o  L 

o 

Y 

=  B,  X’  +  E,  Y  +  H, 

Z* 

* 

L  o  L  o  L 

o 

Z 

* 

=  C,  X’  +  F,  Y'  +  J, 
L  o  L  o  "^L 

Z' 

o 

7.  Applying  ii-3 

=  X.  cos  t  +  Y .  sin  t 
E  * 

Yg  =  -X^  sin  n  t  +  Y*  cos  12  t 
Zg  =  z* 


Or,  we  may  write: 


Xg  =XJAj^  cos  12  t  +  Bg  sin  12  t)  +Y^(Dg  cos  12  t  +  sin  12  t  ) 

+  (Z  +  p)  (G,  cos  12  t  4-  Hj  sin  12  t ) 

O  JU  JLj 

Y„  =X  (- A,  sin  12 1  +  B,  cos  12  t)  +  Y  (-D,  sin  12 1  +  E.  cos  12  t  ) 

Q  Lj  J>^  O 

+  (Z^+  p)  ("CJg  sin  12  t  +  cos  12  t  ) 


=  C,  X  + 
L  o 


F  Y  + 
^  o 


Ir  <Z  + 
^  O 


Pc  ^ 


where: 

A  =  cos  0  cos  (jj  cos  ip  +  sin  0  sin  ip 

Lj  \-i  La  1^ 

Bj  =  sin  0.  cos  <()  cos  tp.  -  cos  0  sin  p 

Cg  =  -sin  cos  p^ 

D,  =  cos  0  cos  <()  sin  p  -  sin  0  cos  - 

E.  =  sin  0  cos  <()  sin  p  +  cos  0  cos  p 

La  Lj  i_(  1j  1j 

F  j  =  -sin  <()  sin  p 

■Lj  .Lj  Lj 

G.  =  cos  0  sin  <() 

■Lj  .Lj  Lj 

Hg  =  sin  0^  sin  <t>j^ 

Jl  =  cos 


4i, 
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A-16 


A-17 


A-18 
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8.  And  applying  ii-4: 

2_„  2  2  2 

^  ^E  ^  ^E 


<1,  =cos-'Zg/yX^2^Yg2+z/ 


A-19 


e  =sin-  +Y^ 

V  Transformation  of  a  Vector  From  the  System  Fixed  in  the  Earth  to  the  Inertial 
System  (x,  y,  z) 

From  Fig.  5,  we  see  that  a  vector  may  be  expressed  as; 


^  ^  ^'eS  ^  ^  ES 


We  wish  to  express  this  vector  in  terms  of  its  components  in  the  i'  „  . 


j'oS  ’  ^'oS  jbH,)  directions. 

This  may  be  done  by  3  successive  transformations; 


1.  From  the  T'^g,  (or  i' ^  ,  k'^,  )  to  the  ,  fg  ,  kg  componaits. 

Applying  ii-1 


'‘'e'’‘e*e-"%®e  +  4'=e 

‘^'e  "  ^E  ^E  ^  ^E  “e  ^‘^E-^E 


A- 20 


/N  'N  A. 

2.  From  the  i„,  j_  ,  k  to  the  i^  j*  k*  components.  Applying  ii- 3: 

£l  £l  £l  < 


A  /N 

i„  =  i  cos  S2 1  +  sin  S2  t 

Ii  ^  ^ 

A  /S  /N 

=  -i  sin  S2  t  +  j  cos  J2  t 

li  ♦  * 

A  A 

^E  ^ 
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^  /N 


*  *  ^  **  T  '  J 

*  o  L  o 

/V 

i  =  i’  D'-  +  j' 
o  L  o 

/\ 


/V  ^ 
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But 


fs.  'y  y\ 

^'oS  "  ^'o  "  W 


to  the  i'  ,  j'  , 
o  •'  O 

k'^  components.  Applying  ii-1: 

B’,  +  k’  C\ 

L.  o  L 

E',  +  k'  F’, 

L  o  L 

A-22 

H',+  k* 

Li  O  ■'  L 

/K  XS 

^  oS  ^  o  ”  ’ 

A  /N 

k'  e  - 

oS  O  0 

Then  the  vector  v  may  be  written: 


V  =Pi'g  +Qk'g 

/V  /N  /X 

=  Ri'  „+  Sj'  _+Tk*  ^ 
oS  •’  oS  oS 

/V  /\  /N 

=  Rlo+SJ^+Tk„ 


A-23 


where: 


R  =  P(Ag  (A'j^  cos  J2  t  +  D’^  sin  n  t)  +  Bg  ("A'^  sin  n  t  +  D’^  cos  Q  t) 

+  Ce  G'l) 

+  Q  (Gg  (A'g  cos  n  t  +  D'g  sin  n  t  )  +  Hg  (-A'g  sin  n  t  +  D'g  cos  n  t) 

+  JeG’l) 

S  =  P  (Ag  (B'g  cos  J2  t  +  E'g  sin  n  t  )  +  Bg  (-B'g  sin  n  t  +  E'g  cos  n  t  ) 

+  CeA’l) 

+  Q  (Gg  (B'g  cos  J2  t  +  E'g  sin  n  t  )  +  Hg  ("B’g  sin  J2  t  +  E'g  cos  J2 1  ) 


+  Je 


A-24 


T  =  P(Ag  (C^  cos  n  t  +  F’j^  sin  n  t )  +  (-C^^  sin  n  t  +  cos  n  t ) 
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and: 


^CeJ'l) 


+  Q  (Gg  (C’j^  cos  J2  t+  F'j^  sin  J2  t )  +  Hg  (-C'j^  sin  n  t  +  F cos  n  t  ) 


Ag  =  cos  0g,  cos  4)g  cos 
Bg  =  sin  0g  cos  <{)jg  cos 
Cg  =  -sin  <{)g  cos  ip^. 


+  sin  0_  sin 
E 

-  cos  0_  sin  tp„ 
E  ^E 


H. 


cos  0„  sin  <{)„ 

Jbi  t 

sin  0g  sin  <{)g 


Je  =‘:os't.E 
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have  been  defined  in  (13  ). 
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